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摘要 :草地 量 (Zoxostege sticticalis L.) 是 我 国 华北 、 东 北 和 西北 地 区 农 牧 业 生 产 的 重大 害虫 ,其 每 年 暴发 成 灾 的 程度 及 造成 的 经 
济 损失 与 越冬 虫 源 基 数 或 面积 密切 相关 ,但 我 国 草地 旦 主要 越冬 区 的 变化 规律 至 今 尚 未 见报 道 c oe 
下 我 国 草 地 螟 越 冬 区 的 变化 规律 ,以 1951 一 2000 年 我 国 草 地 曝 越 冬 场所 及 面积 为 基础 ,结合 政府 间 气 候 变 化 专门 委员 

(IPCC) 第 五 次 评估 报告 发 布 的 BCC-CSM1.1 气候 模式 数据 ,采用 最 大 炉 ( Maxent) 模型 预测 可 4 种 气候 变化 (RCP 2.6, 4.5, 6.0 
和 8.5) 情 景 下 2050s 及 2070s 我 国 草地 量 越 冬 区 的 变化 情况 。 研 究 结果 :1 ) 训 练 数据 售 和 测试 数据 集 的 受 试 者 工作 特征 曲线 
下 的 面积 (AUC) 分 别 为 0.989 和 0.987 ,表明 模型 的 模拟 精度 很 好 。2) 经 Jackknife 方法 检验 ,降雨 量变 化 方差 bio_15 是 最 重要 
的 变量 ,最 湿 月 份 降雨 量 bio_13 次 之 ,最 后 为 10 月 份 降 水 量 prec_10 及 年 温 变 化 范围 bio_7。 这 些 环境 变量 对 模型 的 贡献 率 均 
超过 10% ,是 模型 构建 最 重要 的 环境 变量 。3) 在 各 种 气候 变化 情景 下 ,到 W2050s 及 2070s 我 国 草地 旦 的 越冬 区 面积 和 位 置 相 对 
于 当前 都 有 不 同 程度 的 扩大 和 北 移 。 其 中 高 适宜 越冬 区 面积 为 当前 的 1.41 一 2.94 倍 ,其 质心 位 置 向 北 移动 78.79 一 226.97 km。 

这 些 结果 表明 ,我国 未 来 草地 量 越 冬 场所 将 会 扩大 和 北 移 。 
关键 词 : 草 地 蜡 ;Maxent 模型 ;越冬 区 预测 ;气候 变化 
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Abstract: Thebeet webworm, Loxostege sticticalis L. ( Lepidoptera: Crambidae) is a destructive insect pest of crops and 
fodder plants in northern China. The degree of annual outbreaks of the larval population during the breeding season is highly 
relatedto the size and location of the overwintering area for hibernating larvae, but how this relationship will change in the 
future rémains unknown. Therefore, variations in the size and location of the overwintering area of L. sticticalis were studied 
to reveal changes under various scenarios of global climate change. This was accomplished using overwintering data of L. 
sticticalis in China from 1951 to 2000, combined with the application of the Maxent model and the BCC-CSMI1.1 climate 
change modeling data under four scenarios of RCP , 2.6, 4.5, 6.0 and 8.5, released by IPCC AR5. Our results showed that 


the AUC (area under the receiver operating characteristic curve ) value for the training and test data were 0.989 and 0.987 ， 
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respectively, indicating the data fit the model very well. In addition, the most important variable in the model as determined 
by the jackknife test was the variance in the precipitation coefficient (bio_15 ) , followed by precipitation during the wettest 
month (bio_13 ) precipitation in October (bio_10), and annual temperature range (bio_7). Their individual contribution 
to the model was greater than 10% each. Finally, the centroid of the highly suitable overwintering location would be shifted 
104.93 一 207.47 km and 78.79—226.97 km northward, whereas the size of the area would expanded to 1.41 一 2.00 and 
1.94 一 2.94 fold by the 2050s and 2070s, respectively, compared to those at the present time. 


Key Words: Loxostege sticticalis; Maxent model; overwintering acreage and location forecasting; climate change 


草地 量 是 我 国 华 北 ,东北 和 西北 (三 北 ) 地 区 一 种 重要 的 农 牧 业 害虫 , 曾 多 次 暴发 成 灾 并 造成 重大 经 济 损 
失 ''" 习 。 从 建国 以 来 至 2000 年 我 国 草地 暇 的 主要 越冬 区 分 布 在 38" 一 43*N ,108。 一 118" 玉 之 间 的 农 牧 交错 地 
区 ,这 片区 域 海拔 较 高 (1000 一 1600 m) ,而 从 行政 区 划 方 面 来 看 ,主要 溃 括 山西 北部 铅 北部 和 内 蒙古 中 西 
部 地 区 2] 。 而 在 东北 草地 旦 一 代 常 发 区 , 某 些 年 份 也 有 一 定数 量 的 越冬 虫 源 5 但 油 有 形成 稳定 的 越冬 
工区 !“9 。 草 地 量 必 须 以 老 熟 幼 虫 在 土壤 中 做 莉 滞 育 的 方式 来 越过 严寒 的 冬季 , 这 是 其 年 生活 史 中 较 薄弱 的 一 
c 环 , 因 而 能 否 顺利 越冬 及 越冬 基数 的 大 小 将 会 直接 影响 到 来 年 草地 蜡 是 否 发 生 及 发 生 的 程度 。 罗 礼 智 等 根据 
2003 年 我 国 草地 蜡 2 幼虫 代 发 生 严 重 ,并 在 较 广 阔 的 区 域 留 下 大 量 越冬 喘 源 的 情形 ,成 功 预测 了 2004 年 我 
国 草地 旦 一 代 幼虫 的 大 发 生 !” 。 而 2005 年 我 国 草地 蜡 二 代 幼 虫 发 此 较 辆 面积 较 小 ,因而 该 年 草地 量 的 越冬 
CO 区 越冬 基数 较 小 ,根据 这 种 形势 , 罗 礼 智 等 预测 2005 年 我 国 草地 旦 一 代 幼虫 发 生 程度 较 轻 ,和 当年 实际 的 发 
Ja 生 情况 相符 "中 。 通 过 生产 实践 已 经 证 明 我 国 草地 量 一 代 幼 虫 的 发 生 程度 和 越冬 区 以 及 越冬 基数 的 大 小 有 着 
- ”密切 的 联系 ,因而 能 否 准 确 的 预测 草地 由 越 和 区 和 越冬 磊 数 的 大 四 将 在 一 定 程度 上 决定 能 否 准确 预测 草地 由 
二 一代 幼 虫 的 发 生 程 度 。 而 目前 国内 对 草地 量 越 冬 区 的 预测 却 鲜 有 研究 。 
C 随 着 全 球 气 候 的 变 上 暧 ,我 国 草 地 晨 的 越冬 区 也 有 了 一 - 些 新 的 变化 ,2008 年 草地 量 二 代 幼 虫 在 我 国 三 北 地 
~ 区 大 暴发 ,其 发 生 范围 和 发 生 程度 都 为 历年 历次 滋 最 ,而 随后 二 代 幼 虫 留 下 的 大 量 虫 源 又 成 功 在 东北 地 区 越 
~ 冬 , 并 成 为 次 年 当地 越冬 代 成 虫 大 发 生 的 生 要 电源 。 同 样 是 二 代 幼虫 大 发 生 的 次 年 ,2009 年 春 在 新 疆 和 田 
(5 ”地 区 海拔 4000 m 以 上 的 昆仑 山脉 调 汰 到 永 草 地 蜡 越 冬 幼虫 ""。 这 些 变 化 和 气候 变化 应 该 存在 某 种 联系 。 草 
地 蜡 的 越冬 地 是 草地 蜡 种 群 的 庇护 所 S 基 适宜 的 环境 条 件 为 草地 旦 顺利 度 过 冬季 提供 了 保证 。 而 同时 草地 量 
= 一 的 越冬 地 也 将 是 草地 量 幼 虫 任 长 发 育 繁殖 的 较为 适宜 的 场所 。 越 冬 区 的 大 小 和 位 置 在 一 定 程度 上 会 影响 到 
草地 量 越 冬 虫 源 的 数量 和 里 量 , 欠 而 进一步 影响 的 次 年 草地 晶 发 生 的 数量 和 位 置 。 随 着 气候 的 变化 ,世界 范 
围 内 很 多 物种 都 有 相应 的 变化 ,较为 一 致 的 结论 是 随 着 气候 的 变 暖 ,物种 的 分 布 区 向 着 高 纬度 以 及 高 海拔 的 
方向 移动 玉 轨 。 草 地 嵌 在 我 国 的 越冬 区 呈现 的 这 些 变化 和 其 他 物种 的 变化 趋势 较为 一 致 。 而 随 着 气候 的 进 
一 步 变化 ,我 国 草地 量 的 越冬 区 将 会 如 何 变化 ? 这 对 我 国 未 来 草地 蜡 的 防治 工作 将 会 有 重大 的 影响 。 本 文 将 
结合 目前 通用 的 物种 分 布 模型 ( Maxent) 和 不 同 气候 变化 情景 数据 对 我 国 草地 螟 越冬 区 在 21 世纪 50 年 代 及 
70 年 代 的 变化 情况 进行 研究 。 


0 材料 与 方法 


1 本 草地 量 越 冬 地 分 布 数据 点 来 源 

草地 量 越 冬 地 分 布 数据 收集 于 目前 发 表 过 的 期 刊 论文 ,会 议论 文 ,学 位 论文 ,相关 省 地 植保 站 及 全 国 农 技 
推广 中 心 相 关 草 地 量 越 冬 调查 文件 等 ?45252324 。 在 越冬 地 分 布 数据 处 理 中 去 除了 部 分 范围 较 大 ,模糊 不 清 
的 ,重点 采用 具有 明确 经 纬度 标记 ,以 及 位 置 能 精确 到 乡镇 一 级 的 记录 ,并 通过 这 些 地 名 在 谷歌 地 图 上 查 出 相 
应 位 置 , 并 利用 其 提供 的 经 纬度 工具 查询 记录 地 所 在 位 置 的 经 纬度 。 并 将 经 度 及 纬度 两 者 之 间 相 差 都 不 超过 
0.16666667° 的 记录 点 只 保留 经 度 及 纬度 最 小 的 点 ,以 避免 在 分 析 时 同一 栅 格 内 出 现 多 个 记录 点 。 经 过 整理 
最 终 获 得 114 个 有 效 的 草地 师 主 要 越冬 记录 数据 点 (图 1)。 
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1.2 ”环境 数据 来 源 

本 文 根 据 草地 旦 的 生物 学 习性 ,选取 了 草地 晶 越 冬 
期 间 1 一 4 月 及 9 一 12 月 的 温度 (平均 最 高 气温 ,平均 最 
低 气温 及 平均 气温 ) 和 降水 的 相关 气候 因子 以 及 如 表 1 
所 示 的 19 个 生物 气候 学 因子 .以 及 海拔 高 度 做 为 模型 
预测 的 环境 数据 。1951 一 2000 年 50 年 间 每 月 的 平均 
气温 ,平均 最 高 气温 ,平均 最 低 气 温 以 及 平均 降水 量 , 海 
拔高 度 等 及 19 个 生物 气候 学 指标 作为 草地 量 当 前 越冬 
区 的 环境 数据 。2041 一 2060 和 2061 一 2080 两 个 未 来 
时 间 段 的 相应 气候 指标 作为 预测 越冬 区 分 布 的 环境 数 
据 。 所 有 的 环境 数据 均 下 载 于 Worldclim (http://www. 
worldclim.org) 网 站 ,各 变量 的 空间 分 辩 率 为 10 arc-min 
(=18km) 。2050s 和 2070s 两 个 未 来 时 间 段 的 相应 环 
境 数据 采用 政府 间 气 候 变 化 专门 委员 会 (IPCC ) 第 五 次 
气候 评估 报告 发 布 的 BCC-CSM1.1 气候 模式 4 种 不 同 
气候 变化 情景 (RCP 2.6、RCP 4.5、RCP 6.0 和 RCP 8.5) 
数据 。RCPs ( Representative Concentration Pathways ) 情 
景 为 了 CC 第 5 次 评估 报告 开发 的 以 本 世纪 未 的 辐射 强 
迫 大 小 命名 的 气候 变化 新 情景 , 共 包括 4 种 典型 情景 
(RCP 2.6、RCP 4.5、RCP 6.0 和 RCP 8.5,RCP 2.6 代表 
到 2100 年 辐射 强迫 大 小 达到 2.6W/m 的 温室 气体 排放 


图 1 我国 草 地 紧 主 要 越冬 区 越冬 点 分 布 情况 ( 图 中 的 小 黑 点 代 
表 越 冬 点 》 
Fig. Distribution of overwintering points for L. sticticalis in 


China (dark dots in the figure represent overwintering points) 


量 ,其 他 以 此 类 推 , 分 别 代表 温室 气体 排放 量 低 , 只 低 ,中 高 和 高 4 种 情景 ) 。 与 前 人 研究 常 采用 的 温室 气体 
排放 情景 (Special Report on Emissions Scenaiios ,ISRES) 相 比 ,RCPs 情景 进一步 考虑 了 应 对 气候 变化 的 各 种 策 
略 对 未 来 温室 气体 排放 的 影响 ,因而 更 科学 地 描述 了 未 来 气候 变化 的 预 估 结 


表 1 相关 环境 变量 名 称 及 描述 


Table 1 Description of environmental variables 


变量 Variables 描述 Descriptio 己 变量 Variables 描述 Descriptions 
bio_1 年 平均 气温 bio_2 昼夜 温差 月 平均 值 
bio_ 3 层 夜 温差 与 年 温差 比值 bio_4 温度 变化 方差 
bio_5 最 热 月 份 最 高 温 bio_6 最 冷 月 份 最 低温 
bio 7 年 温 变化 范围 bio_8 最 湿 季 度 平均 温度 
bio_9 最 干 季度 平均 温度 bio_10 最 暖 季度 平均 温度 
bio_ Ml 最 冷 季 度 平均 温度 bio_12 年 平均 降雨 量 
bios13 最 湿 月 份 降雨 量 bio_14 最 干 月 份 降雨 量 
bio_15 降雨 量变 化 方差 bio_16 最 湿 季 度 降 雨量 
bio M7 最 干 季度 降雨 量 bio_18 最 暖 季 度 降 雨量 
bio_19 最 冷 季 度 降 雨量 alt 海拔 

tmax_M M 月 平均 最 高 温度 tmin_M M 月 平均 最 低温 度 
tmean_M M 月 平均 温度 prec_M M 月 平均 降水 量 


2.3 ”Maxent 预测 模型 构建 


Maxent 模型 (最 大 燃 模 型 Maxent 3.3.3 版 本 http://www.cs.princeton.edu/schapire/maxent) 是 Phillips 等 
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开发 用 于 模拟 给 定 环境 条 件 下 物种 出 现 的 概率 的 物种 分 布 模型 '” 。 模 型 输出 为 一 系列 概率 值 的 栅 格 ,概率 
越 大 ,表明 该 栅 格 环境 条 件 越 适 宜 于 物种 生存 。 最 大 炉 原 理 是 根据 已 知 物种 分 布点 对 被 预测 地 域 该 物种 适 生 
性 作出 最 客观 的 推断 ,而 并 不 需要 无 物种 分 布 的 数据 点 。 而 在 历史 的 记录 之 中 ,有 物种 分 布 的 数据 点 能 够 引 
起 人 们 的 高 度 关注 ,而 无 物种 分 布 的 情况 则 相反 ,因此 物种 分 布 的 数据 往往 都 是 只 有 有 物种 分 布 的 存在 点 的 
数据 。Maxent 模型 能 够 用 这 种 只 有 存在 点 的 物种 分 布 数据 来 对 拟 预测 的 区 域 该 物种 存在 概率 进行 很 好 的 预 
测 。 模 型 基于 Java 语言 开发 ,输入 草地 量 越 冬 区 分 布点 文件 ,并 输入 ASCII 格式 的 相关 环境 因子 ,设置 模型 参 
数 , 便 能 运行 模型 ,得 到 相关 参数 。 模 型 还 提供 刀 切 法 (Jackknife) 检 验 , 能 对 环境 因子 贡献 率 和 重要 性 进行 分 
析 , 并 利用 ROC 曲线 ( 受 试 者 工作 特征 曲线 ) 评 价 模型 模拟 结果 的 精准 度 ,ROC 曲线 下 的 面积 (AUC ) 莽 大 代 
表 模 型 模拟 效果 越 好 '*1 。 不 同 的 AUC 取 值 对 应 的 含义 如 下 ;0.5 一 0.6 预测 失败 ,0.6 一 0.7 预测 效果 差 ,0.7 一 
0.8 预测 效果 一 般 ,0.8 一 0.9 预测 效果 较 好 ,0.9 一 1 预测 效果 极 好 。 大 量 实践 证 实 Maxent 模型 是 六 种 用 于 评 
价 物种 栖息 及 越冬 地 适宜 性 效果 较 好 的 物种 分 布 模型 3] 。 

将 草地 蜡 越 冬 区 分 布点 数据 整理 成 包含 经 纬度 的 csv 文件 ,导入 Maxent 模型 ,六 将 1951 一 2000 年 期 间 
bio_1-bio_19 共 19 个 生物 气候 学 因子 1 一 4 月 及 9 一 12 月 的 月 平均 气温 ,月 平均 最 高 气温 ,月 平均 最 低 气温 ， 
月 平均 降水 量 及 海拔 高 度 等 数据 导入 Maxent 模型 的 环境 数据 项 ,再 将 不 同 气候 变化 情景 下 2050s 及 2070s 相 
应 的 数据 导入 Maxent 模型 的 预测 数据 项 ,选中 绘制 环境 变量 响应 曲线 和 环境 变量 重要 性 Jackknife 检验 选项 ， 
并 设置 用 于 模型 随机 检验 的 分 布点 数据 为 20% ,其 它 参 数 保持 默认 值 浊 设置 好 输出 数据 的 相关 路 径 后 运行 
模型 。 

2.4 ”Maxent 模型 预测 效果 评价 及 不 同 气候 变化 情景 下 我 国 草地 蜗 越 冬 区 变化 情况 

对 Maxent 模型 输出 数据 进行 评价 ,根据 训练 样本 及 测试 样本 的 AUC 值 判断 模型 模拟 预测 效果 的 好 坏 。 
考察 各 个 环境 变量 对 模型 的 贡献 率 及 重要 性 。 对 于 输出 的 草地 恤 越 冬 区 适合 性 概率 值 栅 格 图 层 ,选取 特异 性 
和 灵敏 度 之 和 最 大 时 的 概率 值 作为 区 分 草地 螟 越 冬 区 适 宝 坊 否 的 阅 值 (TH, 经 过 计算 本 研究 中 TH= 0.107) 。 
在 ArcGIS 10.1 中 对 所 得 草地 坚 越 冬 区 适合 性 栅 格 图 层 运 用 栅 格 分 类 工具 Reclass 按照 0-TH、TH-0.5 及 0.5 一 
1 分 类 成 非 适宜 区 .适宜 区 及 高 适宜 区 3 类 会 再 运用 分 类 统计 工具 Zonal 计算 3 类 分 区 的 面积 以 及 质心 位 置 。 
比较 适宜 区 及 高 适宜 区 在 不 同年 代 及 不 同 气候 变化 情景 下 面积 和 质心 位 置 变化 情况 。 


3 结果 与 分 析 


3.1 气候 变化 情况 

对 本 文 所 用 未 来 气候 预测 弄 拟 的 BCC-CSM1.1 气候 模式 数据 检测 发 现 , 无 论 是 2050s 还 是 2070s ,我 国 区 
域内 的 气候 相对 于 基准 年 4951 一 2000 年 的 都 有 所 改变 。 随 着 温室 气体 排 量 的 增加 ,2050s 年 均 温 上 升 1.5 一 
3.5% 不 等 ;2070s 奎 均 温 上 升幅 度 在 1.5 一 4% 不 等 ,在 同一 温室 气体 排 量 下 2070s 气温 上 升 的 幅度 较 2050s 要 
大 ;而 北方 地 区 特别 是 东北 地 区 ,年 均 温 的 上 升幅 度 高 于 南方 地 区 (图 2) 。 而 年 降水 量 的 变化 则 较 不 规则 ,全 
国 大 部 分 地 网 在 2050s 及 2070s 年 降水 量 将 会 有 所 增加 ,但 是 新 疆 地 区 ,华北 平原 区 以 及 西南 地 区 则 在 一 定 
程度 正 年 降水 量 将 会 减少 (图 3) 。 
3620 MaXenit 模型 的 预测 评价 

如 图 4 所 示 , 训 练 数据 集 和 测试 数据 集 ROC 曲线 下 的 面积 AUC 分 别 为 0.989 和 0.987, 和 1 非常 接近 ,在 
模型 的 建立 过 程 中 ,随机 选取 了 80% 的 数据 点 作为 训练 数据 集 ,建立 了 预测 模型 ,对 未 参与 模型 构建 的 20% 
的 数据 点 进行 检验 ,发 现 其 符合 度 较 高 ,表明 模型 的 模拟 效果 较 好 。jJackknife 的 检验 结果 表明 ,在 去 除 一 个 变 
量 后 “gain” 值 下 降 最 多 的 是 降雨 量变 化 方差 bio_15 ,也 就 是 说 bio_15 变量 中 包含 了 更 多 其 他 变量 没有 的 信 
息 , 是 最 重要 的 变量 (图 5)。 其 次 ,最 湿 月 份 降雨 量 bio_13 .10 月 份 降水 量 prec_10 及 年 温 变化 范围 bio_7 三 
个 环境 变量 在 模型 中 的 贡献 率 也 较 大 ,都 超过 了 10% , 而 它们 的 变化 情况 对 成 为 草地 量 越 冬 区 的 概率 的 影 
响 呈 现 峰 形 曲线 (图 9) ; 而 10 月 份 平 均 最 低温 tmin_10 温度 变化 方差 bio_4 \ 海 拔高 度 alt 、 最 暧 季度 降雨 量 
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表 2 不 同 环境 变量 在 模型 构建 中 的 贡献 率 与 重要 性 


Table 2 Percent contribution and permutation importance for environmental variables 


环境 变量 贡献 率 重要 性 环境 变量 贡献 率 重要 性 
Variable Percent contribution ”Permutation importance || Variable Percent contribution Permutation importance 
alt 5,7 11.9 prec_4 0.2 1.2 
bio_1 0 0 prec_9 0.1 13.3 
bio_10 0 0.6 tmax_1 0.9 0.3 
bio_11 0 0 tmax_10 0 0 
bio_12 0 0 tmax_11 0 0.2 
i0_13 12.3 17.6 tmax_12 0 0 
10_14 0 0 tmax_2 0 0 
io_15 19.3 14.6 tmax_3 2.6 3.8 
io_16 1.3 0.3 tmax_4 0 0 
bio_17 1.8 13.1 tmax_9 0 0 
bio_18 4.3 0.5 tmean_l1 0.2 0.1 
bio_19 0.3 0.2 tmean_10 0 0 
bio_2 0 0 tmean_11 02 0.1 
10_3 0 0 tmean_12 0 0.1 
i10_4 7.4 0.5 tmean_2 0.2 0.1 
10_5 0.2 0.7 tmean_3 0 0.9 
bio_6 0.1 0.1 tmean_4 0 0 
bio_7 10.1 0.9 tmean_9 0 0 
bio_8 1.1 15.7 tmin_l 0 0 
bio_9 0 0 tmin_10 8.8 0 
prec_1 0 0 tmin_114 0 0 
prec_10 12 0.4 tmins 12 0.1 0.8 
prec_11 0 0.2 tmiaf 2 2.6 0.9 
prec_12 0 0.2 tmiin_3 3.7 0.2 
prec_2 4.2 0 tmin_ 4 0.2 0.3 
prec_3 0 0 tmin_9 0 0.2 


bio_18 二 月 份 降水 量 prec_2 3 月 份 平均 最 低温 tmin_ 
3 3 月 份 平均 最 高 温 tmax_3 .2 月 份 平均 最 低温 tmin_2、 
最 二 季度 降雨 量 bio_17、 最 湿 季 疫 降 雨量 bio_16 和 最 

湿 季度 平均 温度 bio_8, 等 11 从 环境 变量 的 贡献 率 也 都 
超过 了 1% ,这 些 变量 是 模型 构建 比较 重要 的 环境 变量 
( 表 2)。 

刀 切 法 检验 的 步 又 是 先 去 除 某 个 变量 ,再 计算 剩余 
变量 构建 模型 的 (gaim” 值 ,以 及 计算 只 有 这 个 变量 时 构 
建 模型 的 fgain” 值 ,同时 还 计算 包括 所 有 变量 时 构建 模 
型 的 tain 了 值 。 通过 比较 这 些 “gain” 值 的 变化 情况 来 02 
确定 该 变量 的 重要 性 ,去 除 某 个 变量 后 “gain” 值 减少 越 


1.0 


> 
ey 


灵敏 度 Sensitivity 
© 
及 


一 训练 数据 (AUC = 0.989) 
一 测试 数据 (AUC = 0.987) 


多 ,该 变量 在 模型 构建 中 作用 越 大 ,贡献 率 越 高 6 一 随机 数据 (AUC = 0.5) 

3.3 气候 变 化 下 草地 量 越 冬 区 空间 格局 变化 0 02 04 0.6 08 1.0 
Maxent 模型 输出 的 我 国 草 地 量 当 前 越冬 区 位 于 山 Be 

西北 部 河北 张家口 以 及 内 蒙古 乌 兰 察 布 盟 (图 6) ;这 图 4 Maxent 模型 ROC 曲线 

和 我 国 草地 旦 目前 越冬 区 的 范围 是 一 致 的 。 而 在 不 同 Fig.4 ROC curve of Maxent model 

气候 变化 情景 下 ,我 国 草 地 蜡 的 越冬 区 位 置 相 对 于 当前 


都 有 一 定 程度 的 北 移 。 在 2050s RCP 4.5 及 RCP 6.0 情景 模式 下 我 国 草地 临 越冬 区 将 会 完全 移出 山西 及 河北 
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北部 ,并 全 部 落 在 内 蒙古 境内 ;而 2070s 我 国 草地 量 越 9 去 除 某 一 变量 a 只 包含 某 一 变量 ”0 包含 全 部 变量 
冬 区 在 不 同 气候 变化 情景 下 相对 于 2050s 面积 有 所 扩 alt 

大 ,特别 是 高 适宜 越冬 区 的 面积 扩张 明显 (图 6) 。 由 于 op 

越冬 区 的 边缘 并 不 规则 和 整齐 ,因而 采用 越冬 区 的 质心 be- 

来 表征 其 位 置 的 变迁 更 形象 直观 和 具有 代表 性 。 对 草 pt 

地 蜡 的 适宜 越冬 区 质心 变化 分 析 表 明 :2050s 及 2070s bio 6 

我 国 草地 量 适 宜 越冬 区 的 质心 在 未 来 不 同 气候 变化 情 io 

景 下 较 当 前 越冬 区 的 质心 位 置 都 有 一 定 程 度 的 北 移 ,经 区 3 

过 计算 其 经 向 距离 为 93.08 一 230.42 km 不 等 ;其 中 

2070s RCP 2.6 情景 下 草地 蝇 适 宜 越冬 区 质心 经 向 北 移 bi 

距离 最 短 为 93.08 km; 而 2050s RCP 8.5 情景 下 草地 是 Bi-9 

适宜 越冬 区 质心 经 向 北 移 距 离 最 长 为 230.42 km; 同一 pe 

气候 变化 情景 下 草地 蜡 适 宜 越 冬 区 质心 2050s 较 2070s pree 1 

位 置 更 靠 北 (RCP 6.0 除外 ,图 7 B)。 对 草地 蜡 的 高 适 3 Eee 


宜 越 冬 区 质心 变化 分 析 表 明 :不 同 气候 变化 情景 下 我 国 prec 9 
草地 量 高 适宜 越冬 区 质心 在 2050s 及 2070s 相对 于 当 Wax ll 


tmax 11 

前 的 高 适宜 越冬 区 质心 位 置 有 78.79 一 226.97 km 不 等 。 me 
的 经 向 北 移 ;2070s RCP 2.6 情景 下 高 适宜 越冬 区 的 质 ee 
心 相 对 于 当前 位 置 移动 距离 最 小 ,为 78.79 km; 而 2070s we 
RCP 8.5 情景 下 高 适宜 越冬 区 的 质心 相对 于 当前 位 置 。 wall 
移动 距离 最 大 ,为 226.97 km( 图 7 A)。 tmean 2 
3.4 气候 变化 下 草地 量 越 冬 区 面积 变化 on 
如 图 8 所 示 , 随 着 温室 气体 排 量 的 上 逢 ;不 同 气候 tmin 1 


变化 情景 下 我 国 草地 量 的 越冬 区 面积 有 所 增长 ,而 高 适 tmin_ll 
宜 越 冬 区 面积 增长 较为 明显 。 当 前 气候 荣 件 下 我 国 草 Pos 
地 量 越 冬 区 的 面积 通过 Maxent 模型 预测 的 是 36.61 万 nih 9 
km? ;其 中 ,高 适宜 越冬 区 12 万 WRm? ,适宜 越冬 区 24.61 


万 km *。 草 地 临 越冬 区 面积 看 W2050s 不 同 气候 变化 情 08 1012141618 20222426283032 
景 下 变化 不 大 ,但 是 高 适宜 越冬 区 的 面积 变化 则 较为 明 sl dm 


显 , 随 着 温室 气体 排 量 的 增加 ,高 适宜 越冬 区 面积 不 断 ee 
增长 ;2050s 在 REP 8 向 气候 变化 情景 下 高 适宜 越冬 区 
面积 最 类 为 允 .08 历 km ,为 当前 的 2.00 倍 ,在 RCP 2.6 
下 最 小 为 1697 用 km ,为 当前 的 1.41 倍 ( 图 8 A) 。2070s 草地 蜡 越 冬 区 面积 在 不 同 气候 变化 情景 下 变化 较 
天 骤 高 伟 出 现在 RCP 6.0 气候 变化 情景 下 ,达到 了 59.92 万 km2 ,是 当前 气候 条 件 下 的 1.64 倍 ; 高 适宜 越冬 
网 面积 的 最 大 值 也 出 现在 RCP 6.0 气候 变化 情景 下 为 35.28 万 km ,是 当前 的 2.94 倍 ,最 小 值 出 现在 RCP 2.6 
下 为 23.33 万 km ,是 当前 的 1.94 倍 。2070s 其 它 气候 变化 情景 下 高 适宜 越冬 区 的 面积 较 当 前 也 有 比较 明显 
的 增长 (图 8 B) 。 


图 5 环境 变量 刀 切 法 检验 


Fig.5 Jackknife test for environmental variables 


4 ”结论 与 讨论 


本 研究 通过 收集 和 整理 我 国 草 地 蜡 目 前 有 记录 的 越冬 点 ,并 作为 Maxent 模型 预测 的 存在 点 数据 ,结合 
IPCC 第 五 次 评估 报告 发 布 的 BCC-CSM1.1 模式 4 种 代表 性 气候 变化 情景 数据 (RCP 2.6、RCP 4.5、RCP 6.0 及 
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2050 RCP 6.0 ， 2050 RCP 8.5 


2070 RCP 6. 


非 适宜 区 
加 适宜 区 
国 国 高 适宜 区 


图 6 不 同 气候 变化 情景 下 2050s 和 2070s 我 国 草地 蝇 越 冬 区 分 布 情况 图 


Fig.6 Distribution map of Overwintering area of L. sticticalis during 2050s and 2070s under different climate change scenario 


A 226.97 B 
A 207.47 A ?30.42 
42°N 

国 169|88 。 国 163.05 A® 174:75 42oN 

国 169.44 | 173.85- 

@@ 143.19 148 78 132.2 
国 111.23 
@ elo49 ®@ 93.08 


示 夫 0 僵 78.79 41° 


111° 112° 113°E 111° 112° 113° 114°E 


@2050RCP2.6 ®@ 2050RCP 4.5 加 2050 RCP 6.0 息 2050 RCP 8.5 六 当前 
2070RCP26 ®@ 2070RCP4.5 国 2070 RCP 6.0 全 2070 RCP 8.5 


图 7 不 同 气候 恬 化 情景 和 2050s 和 2070s 我 国 草地 暇 高 适宜 越冬 区 ( A) 与 适宜 越冬 区 (B) 质心 分 布 情况 图 ( 图 中 质心 位 置 旁边 的 数字 代 
表 该 点 位 置 果 对 于 当前 质心 点 位 置 的 经 向 距离 ) 


Fig.7 。 Centroids distribution of highly suitable overwintering area (A) and suitable overwintering area (B ) for L. sticticalis in China 


duringi205S0s and 2070s under different climate change scenario ( Numbers beside the mark of centroid represent the longitudinal distances 


between the mark and current centroid ) 


RCP”8.5 ) 对 我 国 草 地 坚 在 2050s 及 2070s 越冬 区 的 分 布 情况 进行 了 预测 。 所 建立 的 预测 模型 预测 的 模拟 较 
好 ,对 目前 草地 坚 越 冬 区 的 预测 结果 和 实际 情况 符合 度 较 高 。 而 预测 模型 中 贡献 率 最 大 的 几 个 因子 则 有 降雨 
量变 化 方差 bio_15 .最 湿 月 份 降 雨量 bio_13,10 月 份 降水 量 prec_10 及 年 温 变 化 范围 bio_7; 从 响应 曲线 上 看 ， 
降雨 量变 化 方差 bio_15 在 100 一 120 之 间 时 ,成 为 草地 暇 越冬 区 的 概率 最 大 ,降雨 量变 化 方差 则 是 代表 一 年 
中 降雨 量 的 各 月 份 的 变异 性 ,变化 越 大 ,其 值 越 大 。 最 湿 月 份 降雨 量 bio_13 在 100 mm 左右 时 ,成 为 草地 量 越 
冬 区 的 概率 值 最 大 。 草 地 量 在 我 国 目前 分 布 于 三 北 地 区 ,这 一 带 属 于 半 和 干旱 半 湿润 地 区 ,一 定量 的 雨水 能 够 
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Fig.8 Overwintering area variation of L. sticticalis in China in 2050s and 2070s under different climate change scenario 
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图 9 重要 环境 因子 响应 曲线 


Fig.9 Response curves of some important environmental variables 


促进 
在 暴雨 过 后 草地 量 的 虫口 数量 急剧 降低 ” 。10 


其 奇 主 植物 的 生长 ,从 而 为 其 生长 发 育 提 伐 


t 食 物 。 然 而 降水 量 太 大 反而 对 其 生长 发 育 不 利 ,有 研究 表明 


的 概率 值 较 大 (图 9) ;10 月 份 草地 旦 老 熟 幼虫 已 


水 量 有 助 于 保持 土壤 的 湿度 ,从 而 有 利于 草地 旦 


草地 旦 越冬 区 的 概率 较 大 ,要 出 现 这 种 温度 变化 范围 ， 


蜡 滞 育 幼虫 的 过 冷却 点 大 概 为 -25.84% La ,冬季 的 气温 低 于 这 
死亡 不 能 安全 越冬 的 危险 。 


月 份 降水 量 prec_10 在 0 一 50mm 之 间 时 ,成 为 草地 量 越 冬 区 
经 进入 清 育 状态 进行 越冬 ,并 且 时 间 还 不 久 , 这 时 适量 的 降 


顺利 越冬 。 年 温 变化 范围 
温度 相对 较 高 而 冬季 


ee 


温度 那么 
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bio_7 在 40 一 50% 范围 


草地 蜡 将 面临 体液 大 量 


内 时 ,成 为 


温度 相对 较 低 。 草 地 
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草地 蜡 是 一 种 迁 飞 性 的 农 牧 业 害虫 ,其 扩散 能 力 较 强 ,东北 华北 和 西北 的 广大 区 域 都 有 它们 的 发 生 和 危害 
甚至 成 灾 。 而 其 生活 史 中 也 存在 薄弱 的 一 环 , 那 就 是 必须 靠 清 育 的 方式 来 度 过 严寒 的 冬季 。 而 其 生理 上 的 理 
化 性 质 决 定 其 越冬 的 地 点 冬季 的 温度 不 能 低 于 其 过 冷却 点 也 就 是 -25.84% ; 同时 温度 也 不 能 太 高 ,否则 清 育 
幼虫 可 能 打破 灌 育 而 变 得 较为 脆弱 ,不 能 安全 度 过 冬季 。 草 地 螟 越冬 区 的 形成 是 各 种 因素 综合 作用 的 一 个 结 
果 。 首 先 要 有 足够 的 越冬 虫 源 ,也 就 是 二 代 幼 虫 要 大 量 发 生 ,并 且 在 环境 条 件 的 作用 下 进入 清 育 状态 ,再 者 当 
地 冬季 的 气温 降水 等 条 件 需要 适宜 草地 量 越 冬 幼虫 的 存活 ;还 有 土壤 理化 性 质 ,土壤 中 的 微生物 以 及 草地 旦 
的 天 敌 等 多 种 因素 都 会 对 草地 临 越冬 产生 影响 ,种 种 因素 相互 影响 ,是 一 个 复杂 的 过 程 。 而 Maxent 模型 下 是 
能 将 各 种 复杂 的 因素 融合 起 来 ,从 中 找 出 一 些 共 性 的 因子 从 而 建立 一 个 较为 完善 的 解释 模型 。 

按照 IPCC 第 五 次 评估 报告 发 布 的 2050s 及 2070s BCC-CSM1.1 模式 的 4 种 气候 变化 情景 数据 进行 预测 ， 
我 国 草地 量 越 冬 区 将 会 较 当 前 在 面积 上 有 所 增长 并 且 位 置 上 有 向 北 移动 的 趋势 。 更 大 的 越冬 区 将 能 承载 更 
大 的 越冬 基数 ,因而 在 未 来 草地 量 的 种 群 数量 可 能 较 当 前 有 更 大 增长 ,其 暴发 成 灾 的 可 能 性 会 更 大 ,其 防 控 任 
务 也 将 会 更 加 艰巨 。 
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